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Low energy proton capture reactions have been studied in statistical model using the semi- 
microscopic optical potential in the mass range A=110-125. Nuclear density obtained from rela¬ 
tivistic mean field calculation has been folded with nucleon-nucleon interaction to obtain the optical 
potential. Theoretical results have been compared with experimental measurements to normalize the 
potential. Results have also been compared with a standard calculation available in the literature. 
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Creation of nuclei heavier than iron occurs mainly 
through three processes - slow and rapid neutron cap¬ 
ture processes and p-process Though some studies 

have shown that p-process amounts to only 0 . 1 % of the 
total abundance of nuclei heavier than iron, there are 
possibly 35 isotopes on the neutron deficient side of val¬ 
ley of stability expected to be produced either entirely 
or partly through the p-process. The p-process depends 
on reactions such as (p, 7 ), ( 7 ,p) and ( 7 , 71 ) reactions, as 
well as on the seed nuclei. 

In spite of the high Coulomb barrier, various astro- 
physical sites may produce medium mass proton rich 
nuclei through radiative proton capture reaction with a 
huge proton flux. One of them is the x-ray burst where 
the peak temperature reaches 1-3 GK. In heavier nuclei, 
it is important to study the ( 7 ,p) reactions which destroy 
proton rich nuclei and compete with the ( 7 , n) reactions. 
In terrestrial laboratories, it is easier to study (p, 7 ) re¬ 
actions and infer the cross-sections for the inverse ( 7 ,p) 
reactions. Besides, in astrophysical environments where 
photodisintegration occurs, the temperature is very high 
and nuclei are likely to exist in excited states as distinct 
from our laboratory experiments. Hence, the contribu¬ 
tion of the ground state may not be substantial. All these 
arguments have led to the study of the forward (p, 7 ) re¬ 
actions to constrain the rates. 

To study the actual abundance of different nuclei and 
evolution of the process, a network calculation involving 
a large number of nuclear reactions is required. The net¬ 
work calculation in an explosive astrophysical situation 
has to take many quantities such as temperature, pres¬ 
sure, proton mass fraction, forward and inverse reaction 
rates etc in to account. Thus, we need to tune the in¬ 
teraction potential properly. As the p-process proceeds 
along proton rich side of the stability valley, it involves 
many nuclei which are unstable and inaccessible as tar¬ 
gets on earth. Hence, theory remains the principal guide 
to gather information about various reactions. It is use- 
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ful to construct the nucleon-nucleus potential through 
folding of the NN interaction with the nuclear density. 
Nuclear density for stable targets may be experimentally 
measured by electron scattering. In unstable nuclei ex¬ 
perimental nuclear densities are not available and theo¬ 
retical results may be used. 

Rauscher et al. [ 1 , @ have calculated the reaction rates 
for various proton, neutron and alpha induced reactions 
and their inverse reactions in Hauser-Feshbach formalism 
for various targets over wide ranges of atomic numbers, 
masses and temperatures in global approach. Optical 
potential, a key ingredient for Hauser-Feshbach calcula¬ 
tion, is often taken in a local or global approach. In the 
present work we concentrate on the nuclei in the mass re¬ 
gion A=110-125. There are a number of p-nuclei in and 
around this region, viz. ^^^In, 

I24,i26^g^ Away from the stability valley, reaction rates 
calculated even from semi-microscopic optical potential 
and phenomenological density prescription, become very 
uncertain and have to be varied by a large amount Q . We 
have used a more microscopic approach where an appro¬ 
priate microscopic NN interaction has been folded with 
nuclear density from mean field calculations to reduce 
this uncertainty. This process can be extended to the 
unstable nuclei. It is imperative to test the theoretical 
method in reactions where experimental data are avail¬ 
able to verify its applicability before extending it to un¬ 
known isotopes. Hence, reaction cross-section for avail¬ 
able stable targets in the mass region have been studied 
in the present work. 

The density dependent M3Y interaction (DDM3Y), 
which is known to give satisfactory results in many cases, 
has been chosen for our purpose. Similar studies in mass 
region A=55-100 [Ml have already been carried out 
and our aim is to extend this approach to A=110-125 
region. 

Energy being low in astrophysical environment, the nu¬ 
clear skin plays a very important role in the nuclear re¬ 
actions. In astrophysical situation nuclear reaction rate 
is sensitively dependent on nuclear density profile. Reli¬ 
able density information can be obtained from relativistic 
mean field (RMF) theory, a very useful tool for explaining 
various low energy nuclear pro pert ies and nuclear struc¬ 
ture as well as nuclear density [Ij, [ij . 
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Nuclear density profiles have been calculated using the 
Lagrangian density FSU Gold [T^. This Lagrangian den¬ 
sity contains, apart from the usual terms for nucleons 
and mesons, two additional non-linear meson-meson in¬ 
teraction terms, whose main role is to soften the EOS of 
symmetric matter and reduce the symmetry energy. 

Calculation of density has been carried out in co¬ 
ordinate space assuming spherical symmetry. A pairing 
force of 300 MeV-fm for both protons and neutrons has 
been used under BCS approximation. More details for 
the method are available in Ref. [l^. Considering the 
finite size of the proton, charge density is obtained by 
convoluting the point proton density. 

Pch{r) = e p{vi)g{v - r/)dr/ (1) 


Here, g(r) is the Gaussian form factor given by, 

gif) = {afTT)~^expi-f/f) (2) 

where a is a constant whose value is assigned to be 0.8 
fm. 

It is well known that RMF can provide an excellent de¬ 
scription of nuclear ground state properties in this mass 
region.To check the validity of nuclear density calcula¬ 
tion, nuclear charge radii may be compared with mea¬ 
sured values. In table U calculated binding energies and 
charge radii values have been compared with experimen¬ 
tal measurements. The theoretical binding energy values 
have been corrected following the phenomenological pro¬ 
cedure in Refs. [13 111- Here we have listed all the stable 
isotopes which occur in the concerned mass region and 
whose experimental charge radii values are known. 

The relative difference between theoretical charge radii 
and experimental values is less than 0.5%. Binding en¬ 
ergy values show a maximum error of 0.2%. Experimen¬ 
tal binding energy and charge radii have been taken from 
Audi et al. [l^ and Angeli [l^, respectively. Effect of 
centre of mass correction would be small for heavy nuclei 
as relative correction in radius goes as A“^/^, where A is 
the mass number of the nucleus 211. 

The M3Y interaction [2^ , even in absence of ex¬ 
plicit density and energy dependence, is able to produce 
reasonably satisfactory results for processes like elastic 
scattering and reaction. After incorporation of den¬ 
sity dependence, it is termed as density dependent M3Y 
(DDM3Y) interaction . The effective nucleon-nucleon 
interaction incorporates the density term as 


giP, e) = Cil- Pie)p^/^) 


(3) 


p and e are the nucleonic density and energy per nucleon, 
respectively. The constants C and (3 have been assigned 
to be 2.07 and 1.624 fm^, respectively 1^ obtained from 
optimum fit for a radioactivity data [26l| . Here we also 
have used Scheerbaum spin orbit potential term with 
phenomenological complex potential depths. The values 
used here for these ^enomenological potential depths 
are same as in Refs. [§-[lll|. 


The DDM3Y interaction has been incorporated in 
TALYS1.4 code [2^ for Hauser-Feshbach statistical 
model calculations. We have used the Goriely HFB- 
Skyrme calculations for theoretical nuclear masses in 
TALYS code for the nuclei whose experimental masses 
are not available. The prescriptions for level densities and 
gamma strength function are taken as Constant Temper- 
ature Fermi Gas model and Goriely’s Hybrid model [^ . 
respectively. All these features are available in the code. 
We have included the effect of the width fluctuation cor¬ 
rection which has a significant impact at low incident 
energies. We have considered a maximum of 30 discrete 
levels in target nuclei and in residual nuclides. Full j — I 
coupling has been used. 

We have compared theoretical and experimental as- 
trophysical S-factors. S-factor is a more slowly varying 
function of energy than the cross-section, and useful for 
comparison of results at very low energy. 

The Gamow window for proton capture reaction is cen¬ 
tered around the Gamow energy in MeV 

Eo = 0.1220(^2 (4) 

with width in MeV 

A = 0.2368(Z2z2/iT|)i/® (5) 

Here, Zp and Zt are atomic numbers of the projectile 
and the target, p is the reduced mass in a.m.u. of the 
system and Tg is the temperature in GK. In the reactions 
considered in the present work Gamow window has a 
range from 1.8 to 4.5 MeV for temperature around 2-3 
GK. Experimental masses have been taken from Ref. [l^ 
in converting obtained reaction cross-section values to S- 
factor for comparison with experimental S-factor values. 
For f,p) reaction, it is worth noting that effective energy 
window will be centered at 

E^^f = Eo + Qp^ ( 6 ) 

because photodisintegration requires energy in excess of 
the threshold value, keeping the width unchanged. 

In astrophysical environments temperature may be 
high and nuclei can exist in ground state as well as vari¬ 
ous excited states. Hence, Maxwellian averaging has been 
performed for reaction rates. 

The optical potential has been normalized for the mass 
region considered to fit the experimental observations. 
The DDM3Y is a purely real interaction. We have cho¬ 
sen the real part of the potential by multiplying the fold¬ 
ing potential with a factor 0.9. The imaginary part of 
the potential has been chosen identical with the real one. 
As we will show, it gives a reasonable agreement with 
the proton capture reactions in this mass region. Dif¬ 
ferent normalizations for individual reactions may fit the 
observed data better for various reactions in this mass 
region. However, such a procedure will lead to a de¬ 
crease in predictive power. In extending the calculation 
to unknown nuclei in the mass region, a single set of pa¬ 
rameters will serve the purpose. 
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TABLE I: Calculated binding energy (MeV) and charge radii (fm) of nuclei in A=110-125 region compared with measured values. 


Nucleus Binding energy Charge radius 
Theory Expt. Theory Expt. 

Nucleus 

Binding energy 
Theory Expt. 

Charge radius 
Theory Expt. 

Nucleus 

Binding energy 
Theory Expt. 

Charge radius 
Theory Expt. 

llOpd 

939.78 

940.21 

4.552 

4.578 

“^Sn 

952.46 

953.53 

4.609 

4.594 

i2"Sn 

1047.73 

1049.96 

4.695 

4.676 

i“Cd 

939.24 

940.64 

4.575 

4.574 

“"Sn 

970.64 

971.57 

4.622 

4.610 

i2iSb 

1026.62 

1026.32 

4.690 

4.680 

“^Cd 

947.98 

947.62 

4.581 

4.580 

“"Sn 

979.14 

979.12 

4.630 

4.617 

i23sb 

1042.30 

1042.10 

4.704 

4.688 


956.38 

957.01 

4.588 

4.595 

“®Sn 

987.24 

988.68 

4.638 

4.627 

i22Te 

1033.75 

1034.33 

4.712 

4.708 

“"Cd 

964.06 

963.55 

4.560 

4.601 


995.30 

995.62 

4.645 

4.632 

i23Te 

1041.99 

1041.26 

4.719 

4.711 

ii^Cd 

971.82 

972.60 

4.604 

4.614 

“»Sn 

1003.27 

1004.95 

4.653 

4.641 

i24Te 

1050.13 

1050.69 

4.725 

4.718 

ii«Cd 

986.47 

987.43 

4.619 

4.628 

i«Sn 

1010.96 

1011.43 

4.660 

4.645 

i25Te 

1058.16 

1057.26 

4.732 

4.720 


963.88 

963.09 

4.605 

4.602 

i^°Sn 

1018.57 

1020.54 

4.667 

4.654 

i2"Xe 

1046.43 

1046.26 

4.755 

4.762 

ii®In 

980.06 

979.40 

4.621 

4.617 

Sn 

1033.52 

1035.52 

4.681 

4.666 
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Comparison of theoretical astrophysical 


S-factors with experimental data for the indicated proton capture reactions. 


The comparison of S-factors for a number of reactions 
in this mass region are presented in Fig. [I] We have 
included all experimental work where reasonable amount 
of data exist. In all the figures solid lines denote the theo¬ 
retical results. Open and solid circles with error bars de¬ 
note experimental data. Experimental data for ^^^Sb and 
^^^Sb have been taken from Harrisopulos et al. [^. In 
that work experimental cross-sections were obtained from 
7 -angular distribution measurement with large volume 
HPGe detectors. Experimental data are from Chloupek 
et al. [13 for ^^^Sn and ^^®Sn, Famiano et al. [13 for 
^^"‘Sn and ^^®Sn, and Ozkan et al. [13| for ^^®Sn, respec¬ 
tively. Experimental data for ^^^Sn and ^^®Sn were ob¬ 
tained using activation technique with BGO detector and 
statistical model calculation was used to describe the re¬ 
sults. For ^^"^Sn and ^^®Sn activation technique was used 
to extract the cross-section [13 ■ 

One can see that the present method can describe the 
experimental results except in a few instances. Of course 
the errors are large in some cases. Except for the proton 
capture by ^^®Sn, all the other reactions are reasonably 
described by the present procedure. It is worthwhile to 
note that statistical model calculations by the authors 


failed to reproduce the cross sections for ^^®Sn [s^ - 

We have compared the astrophysical rates calculated 
using NON-SMOKER formalism [1, [13 and our 

present calculation for reactions which involve p-nuclei 
as target or product to see how the present reaction rates 
differ from the standard ones available in the literature. 
In Fig. [21 we plot some of the rates where the differences 
between the two calculations are significant. The only 
exception is the ^^^In(p, 7 ) reaction, which have been 
shown to emphasize that in some cases, the two methods 
agree to a great extent. The temperature Tg is in GK. 
The data have been plotted in the temperature range 
1.5-4 GK relevant to the above said Gamow window. 
We have used the NON-SMOKER results for astrophys¬ 
ical reaction rate data with Finite Range Droplet Model 
(FRDM) [13 nuclear masses for comparison. We find 
that our results are larger than the NON-SMOKER rates 
by a factor of 2-3 approximately in most cases. In some 
other reactions, the two results are very close to each 
other. It should be interesting to study the effect of the 
present rates on the p-process. 

To summarize, a microscopic study of low energy (p, 7 ) 
reactions has been undertaken in the mass range A=110- 
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FIG. 2: Comparison of proton capture rates calculated from NON-SMOKER (dotted line) calculations and present work (solid 
line) for certain reactions involving p-nuclei. 


125. RMF calculation has been performed to obtain the 
nuclear density profile. The real and imaginary parts 
of the optical model potential have been obtained by 
suitably normalizing the folding potential arising out of 
the DDM3Y interaction to explain the observed reaction 
data. S-factor values have been compared with exper¬ 
imental data. Proton capture reaction rates obtained 


from our calculation have been compared with NON- 
SMOKER reaction rates. 
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